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Our goal is to study statistical properies of "dielectric resonances" which are poles 

of conductance of a large random LC network. Such poles are a particular example 



0^ ■ of eigenvalues A„ of matrix pencils H — AW, with W being a positive definite matrix 



and H a random real symmetric one. We first consider spectra of the matrix pencils 
Q_) ■ with independent, identically distributed entries of H. Then we concentrate on an 



infinite-range ("full-connectivity") version of random LC network. In all cases we 
calculate the mean eigenvalue density and the two-point correlation function in the 
Cn ■ framework of Efetov's supersymmetry approach. Fluctuations in spectra turn out 

1/^ ' to be the same as those provided by the Wigner-Dyson theory of usual random 

g : matrices. 
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I. INTRODUCTION 



Large random reactance networks (that is networks made of random mixture of capacitances C 

and inductances L) possess a peculiar property noticed first by Dykhne W: they have a finite real 

r"| ■ conductance and thus can disperse an electric power. Explanation of this somewhat paradoxical 

v^ ' property is simple, however. Indeed, if such a network is large enough there always exist circuits of 

^ [ "resonance type", with (purely imaginary) conductance showing poles at some frequencies. Then 

^> , the real part of the conductance as a function of frequency consists of a set of (5— like peaks at 

^ those resonance frequencies. When the volume of the network grows to infinity, this set becomes 

more and more dense. Then, adding an arbitrary small (infinitesemal, but fixed) active part to all 

inductances ( one may think, e.g. of the inductance L on each bond being in series with a weak 

resistance R) results in a finite active resistance of the network. 

A random mixture of two active conductances is very well studied in correspondence with the 
bond percolation problem, see e.g. [p[- At the same time, the random reactance networks are 
relatively less studied. 

It is necessary to mention that the random LC (more generally, RL — C) networks emerge in 
various physical contexts. As was shown long ago by Shender H the conductance of the random 
LC network turns out to be intimately related to properties of collective excitations in spin glasses. 
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The mapping between the two problems is possible by an analogy noticed first by Kirkpatrick 
between the Kirchhoff 's law and the equation of motion for the spin operators. 

More recently, RL — C networks were claimed to be an adequate model for describing the 
optical absorption in disordered metal films, see B and references therein, which showed some 
unusual features. This fact motivated Luck and collaborators for a series of insightful numerical 
investigations of two-dimensional RL — C arrays P,Q . 

Since it is the conductance poles (resonances) of the LC networks which dominate properties 
of the weakly dissipative RL — C networks, it is natural to try to understand their properties in 
a greater detail. To access those poles it is convenient, following |0,|| to start with the Kirchhoff 
equations for the electric potential at vertices i of a network: 

Y.^^j{V,~V,)^0 (1) 

where the summation goes over all vertices j which are neighbours of a given i. The equations 
Eqs.(|]) should be complimented with boundary conditions at electrodes. For a two-terminal ge- 
ometry one assumes Vi = (resp. Vi — V) for the vertices i belonging to the left (resp. right) 
terminal. 

It is useful to introduce the (positive semidefinite) Laplace operator D on the network by 

(Dl/), = ^(VS - V.) (2) 

j 

with a convention Vi = on both electrodes. 

In a random LC network each conductance at frequency / = uj/2t: is equal to either ctq = iCco 

or (Ti = (iLu))^^, with a specified probabilities (in what follows we concentrate on the case of equal 

probability for finding L and C bonds in the network). Then the Laplace operator can be written 

as a sum Dc +Dl of its components on the L— and C— bond sets, respectively. It is easy to show 

that the poles of the conductance occur at frequencies / given by the roots of the equation MM : 

det (Di - AD) = (3) 

where we introduce the ratio 



<7o~ai 1 + (w/wo)^ 



^=T^^- .r;;'.,. ; o<A<i (4) 



with 0^0 = l/i^LCY^"^ being a characteristic resonant frequency of the LC network. 

The simplest, but very informative way to understand properties of the random binary mixtures 
on average is to write an equation for the mean conductance S using an effective-medium approxi- 
mation (EMA). For a network with a coordination number z and equal concentrations of ctq— and 
(Ti— bonds it reads as follows [Q: 

S — ctq E — fji 



2cro + {z~ 2)11 2(71 + (z - 2)S 



= (5) 



Analysing the solution of this quadratic equation as a function of A in the interval < A < 1, 
one finds that generically for z > 4 the real part of S(A) is non-zero only inside some interval 
^min{z) < A < Xmax{z)- According to the discussion above it means that a mean density of 
resonances is finite only inside that interval. 

Useful and simple as it is, EA^IA suffers from an essential drawback: it systematically neglects 
fluctuations, whereas taking the fluctuations into account could be important. For example, as was 
noted in y,|7|] the existence of sharp edges Xmin,max was an artifact of the EMA approximation. In 
fact, the density of the resonances is never exactly zero inside the whole interval A S [0,1] due to 
the so called Lifshitz tails, which is a purely fluctuation phenomenon. Apart from that, it seems 
hard to construct EMA as a systematic approximation with respect to some small parameter, 
though there are indications that it becomes progressively exact for higher spatial dimensions d M 
and can be an extremely well- working one for many realistic systems [Ej . Therefore, it is not clear 
how to correct it systematically or to include fluctuation effects into account. 

These unsatisfactory features should be contrasted with the status of the mean-field approxima- 
tion in the theory of magnetism, with which EMA shares its main ingredients. As is well known, 
the mean-field equations become exact for the model with the infinite range of spin interaction. 
Even for strongly disordered spin systems like spin glasses the latter model provides an adequate 
basis for constructing a mean-field theory with many non-trivial properties ||lO|| . 

It is natural to try to consider a model of similar type also for the disordered reactance networks. 
This just amounts to considering the Kirchhoff equations Eqs.(|lj) on a disordered full- connectivity 
graph of N vertices (nodes) connected by N{N — l)/2 edges (bonds), each independently taking a 
value (To or ci. 

Actually, we find it convenient to consider a graph with iV+2 nodes, among them two "terminal" 
nodes (labeled A and B, respectively) are singled out by being attached to external voltage, so 
that the potential at A is equal to V, whereas the potential at B is kept zero. The rest of N 
"internal" nodes are labeled by index i = 1, ..., A^ and the corresponding (induced) potentials are 
denoted with Vi. The nodes are connected in a full-connectivity graph of {N + 2){N + l)/2 — 1 
bonds, with a single bond being excluded by obvious reasons: that connecting terminals A and B 
directly, see Fig.l. 



B 



V 



FIG. 1. Schematic view of the fuU-connectivity graph with three internal nodes. One direct bond 
connecting the terminal nodes A and B is excluded. 



We then attribute conductances a^^i, = <J^fj_, pL,v — A, B,l, ..., N to each bond. It is conve- 
nient to introduce three N— component vectors: v = {Vi, ...,Vn)'^ , Q_a — (cai, ..., cr^w)"^, g_g = 
(ctbi; ..., ctbat)^ (we used here T to indicate the transposition) and N x N matrix S with the 
following structure: 



/ 






(TAl + (TBI -r A^k^i 
CT21 



-0'12 



'TA2 + CrB2 



2^k^2 ^2fe 



- 172 AT 



\ 



\ —CTni —CrN2 ■■■ CFAN + <TBN + J2k^N ^Nk J 

In these notations iV + 1 Kirchhof's equations Eq.([|) (N for the internal nodes i and one for the 
node B) can be written, respectively: 



Su = VoiA 



O-rW = / 



(6) 



where / stands for the total current flowing outwards through the terminal node B. This system 
of equations can be readily solved yielding the expression for the network conductance: 



YAB-y^gigi: g_j^ 



(7) 



If we now deal with a binary networks when each aij is cto with a probability p and ui with the 
probability 1 — p it is convenient to introduce the "symmetric" variables hij such that hij — —1 
if Gij ~ ctq and hij = 1 if Oij = cri, so that ct^ = ^ ([co + o"!] + [ci ^ croj/iy ). In terms of these 
variables the network conductance Eq.(0) can be written as: 

1 



Vab = Yab/(To 



1 r 

2A 



\e' 



HAW 



ILa ^ -^fi 



(8) 



where A = 2A — 1, 
following elements: 



(1, 1, ..., 1), h^ — (/i_Bi, ■■■yhsN) and N x N matrices W,H have the 



W^J = (iV + 2)(5y - 1 , H^j = S^J \ hM + h,B+Y^ Kk I - (1 - 5^3)h^j (9) 

With these expressions in hand we see that resonances of our network are determined by values of 
A satisfying the condition: 

det (h - Aw) = (10) 

Of course, it is evident that the fuU-connectivity construction would be completely inadequate 
one for describing the percolation problem which dominates properties of binary mixtures of two 
real conductances. We are, however, interested in studying the resonances of disordered LC net- 
works, and the existence of such resonances is in no way precluded by "all-to-all" geometry. We 
will find, that such a model turns out to be, in essence, exactly soluable in the limit A^ ^ 1. 
We start with deriving the mean density of resonances which turns out to be a function of the 
scaled variable r = XN^^"^. One can envisage such a scaling already from the EMA result Eq.(^, 
which for z 3> 1 gives \Xmin.max — 1/2| oc Xmin.max ^ z~^". At the Same time, in contrast to 
EMA the support of the spectrum in the infinite-range model does not have artificial sharp edges 
Xmin,max, but rathcr the mean spectral density p{r) smoothly decays to zero as long as r -^- cxi as 
p{r) oc e"*" /^. 

More interesting is the fact, that the infinite-range model opens a possibility to study fluctu- 
ations of various quantities. In the present paper we concentrate on spectral fluctuations and, 
correspondingly, study the two-point correlation function of the resonance densities. The latter 
turns out to be essentially the same as given by the famous Wigner-Dyson theory of random matrix 
spectra. This fact favourably agrees with overall numerical results for the resonance spectra of 
two-dimensional disordered LC networks. An origin of relatively small, but noticable deviations 
from the Wigner-Dyson statistics detected in ||7| remains unclear for us at the moment and could 
be related to 2D features of the networks studied there which are not captured adequately by our 
inflnite-ranged model. This issue deserves further investigations. 

Actually, flnding the set of A's satisflng Eq.(||) or Eq.(|lO|) is an example of the generalized 
eigenvalue problem. The combination of matrices D^ — AD or H — AW in this respect is known 
in the mathematical literature as a pencil of matrices [|ll| or just the matrix pencil. Theory of the 
matrix pencils has many important applications such as e.g. vibration and bifurcation analysis in 
complicated structures [Of and game theory |13]. 

At the same time, it seems that the present knowledge on statistical properties of generalized 
eigenvalues of the pencils of random matrices is rather scarce. The paper p4| considers the mean 
number and the density of real eigenvalues for a pencil H — AW, with both H and W being 
matrices with all independent real entries and no symmetry conditions imposed. At the same 
time, our original physical problem motivates us to be interested in the pencils formed by real 
symmetric matrices H, W, with W being positive definite. It is clear that W > ensures all the 
eigenvalues of the matrix pencils to be real (such pencils are sometimes called in the literature the 



regular ones |ll|). Indeed, in that case the generalized eigenproblem: Hx = AWx is equivalent 
to a usual one: W-^l'^WW-^/'^tj_ = Xy , y^ W^^'^x, with H = W^i/^HW-i/^ being a real 
symmetric matrix. 

The mean eigenvalue density for matrices of similar types (when both W and H are random) was 
studied in some generality starting from the work by Marchenko and Pastur, see p5[ . In fact, the 
mean eigenvalue density can be found for the ensemble H, with H being a real symmetric matrix 
with independent, identically distributed (i.i.d.) entries using a generalized version of results by 
Pastur |16| . We are however not aware of any systematic study of spectral correlations of regular 
pencils of the random matrices. 

This should be contrasted with a very intensive research on eigenvalues of the random matrices 
(which is a particular case of W = 1) performed in recent years in the domain of theoretical 
and mathematical physics p7[ . Below we summarize the most important facts known from these 
studies. 

As is well established |lq- |20| , p2| , p3| , the statistical properties of real eigenvalues Xi of large NxN 
self-adjoint random matrices H are to large extent universal, i.e. independent of the details of the 
distributions 7^(11) of their entries. 

It is important to mention the existence, in general, of two different characteristic scales in 
the random matrix spectra: the global one and the local one. The global scale is that on which 
the eigenvalue density, defined as p{X) — -^Tr(5(X — H), changes appreciably with its argument 
X when averaged over 7^(11). For matrices whose spectrum has a finite support in the interval 
X e [A, B] the global scale is just the length of this interval. 

In contrast, the local scale is that determined by the typical separation A = (Xi—Xi-i) between 
neighbouring eigenvalues situated around a point X, with the brackets standing for the statistical 
averaging. It is given therefore by A = {{Np{X)))^^ . If we are interested in those values of X 
that are sufficiently far from the edges of the spectra the global scale is , roughly speaking, by a 
factor of N larger than the local one. In other words, the mean density {p{X)) can be considered 
as a constant one on the scale A. 

The degree of universality is essentially dependent on the chosen scale. 

As to the global scale universality, one can mention first of all that for the matrices with i.i.d. 
entries under quite general conditions (see, e.g., ||24| and references therein) the mean density is 
given in the limit A^ ^ oo by the so-called " Wigner semicircle law" : 



(pW> = ^v/I^^^^^-^ (11) 

In this expression the parameter a just sets the global scale in the sense defined above. It is 
determined by the expectation value a^ — {^TrlP). It is generally accepted to scale the entries 
Hij ~ 1/A^^/^, i.e. in such a way that a stays finite when A^ — > oo, the local spacing between 
eigenvalues in the neighbourhood of the point X being therefore A ex 1/A^. 

From the point of view of universality the semicircular eigenvalue density is not extremely robust. 



Most easily one violates it by considering an important class of the so-called "invariant ensembles" 
characterized by a probability density of the form 7^(H) ex exp— A^TrT^(H), with T^(H) being an 
even polynomial. The corresponding eigenvalue density turns out to be highly nonuniversal and is 
determined by the particular form of the potential V{H). Only for 1^(H) = H^ it is given by the 
semicircular law, Eq.(0). Actually, any "deformation" Hi = Ho + H of the ensemble H with i.i.d. 
entries by a given fixed matrix Hq results in the mean eigenvalue density belonging to a family of 
"deformed semicircular laws" discovered by Pastur jl^ . 

Moreover, one can easily have a non-semicircular eigenvalue density even for real symmetric 
matrices S; Sij = Sji with i.i.d. entries, if one keeps the mean number of non-zero entries p per 
column to be of the order of unity when performing the limit N —> oo. This is a characteristic 
feature of the so-called sparse random matrices |23,p3,pffl characterized by the following probability 
density of a given entry 5"^ : 

r{S,,) = (1 - ^)6iS..,) + ^h{S.,) (12) 

where h{s) — h{—s) is an arbitrary even distribution function satisfying the conditions: h(Q) < 
oo; J h(s)s^ds < oo. 

Remarkably, a much more profound universality emerges for the two-point correlation spectral 
function defined as: 

{p{Xi)piX2))^ = (p(Xi)) 5{Xi - X2) - y2iX,,X2), (13) 

where we defined the connected part of the correlation function in a usual manner: {AB)c = (AB) — 
{A){B). The nontrivial part of the spectral correlator is called the cluster function y2{Xi, X2)- It 
is one of the most informative statistical measures of the spectra, ||28| . It turns out, that already 
the global scale behaviour of y2{Xi,X2) (i.e. one for the distance S being comparable with the 
support of the spectrum in the limit N — > 00) is rather universal. It is the same for all the 
"invariant ensembles" p9| ] and for those with i.i.d. matrix elements p3| and is determined only by 
the positions of edge points [A, B] of the spectrum pTf . 

Even more interesting is the fact that universality of the correlation function Eq. ( p^ ) (as well 
as of all higher correlation functions) extends to the local scale, i.e. for the distances \Xi — X2I 
comparable with A. This fact was rigorously proved for the unitary invariant ensembles [|l8| , 
extended to the unitary "deformed" ensembles |^1|] and heuristically verified for other invariant 
ensembles pG] as well as for the ensembles of sparse matrices [ p2p9| ] . The particular form of the 
cluster function is different from that typical for the global scale. It is dictated in general by global 
symmetries of the random matrices, e.g. if they are complex Hermitian or real symmetric pq |. All 
specific (nonuniversal) properties are encoded in the value of local spacing A. For the distances 
S such that S ^ A local expressions match with the global one, the latter taken at distances 
S'< [A-B]. 

It turns out, that it is the local scale universality that is mostly relevant for real physical systems 
pO{. Namely, the statistics of highly excited bound states of closed quantum chaotic systems of 



quite different microscopic nature turn out to be independent of the microscopic details when sam- 
pled on the energy intervals large in comparison with the mean level separation, but smaller than 
the energy scale related by the Heisenberg uncertainty principle to the relaxation time necessary 
for a classically chaotic system to reach equilibrium in phase space pl| . Moreover, the spectral 
correlation functions turn out to be exactly those which are provided by the theory of large random 
matrices on the local scale p2,R3], with different symmetry classes corresponding to presence or 
absence of the time-reversal symmetry. 

Motivated by lack of general theoretical results for spectral properties of the pencils of random 
matrices, we start in Sec. II with considering an abstract pencil , with real symmetric H belonging 
to the Gaussian ensemble of random matrices with i.i.d. entries and W being a real symmetric 
positive definite one with fixed given entries. We derive an expression for the mean spectral density 
of such a pencil and demonstrate that the correlation properties on the local scale (in the sense 
defined above) are the same as given by the Wigner-Dyson expressions. Then, in Sec. Ill we extend 
our consideration to the case of resonance statistics in the random infinite-range LC-network, by 
deriving the mean density and the two-point spectral correlation function for the resonances. 

II. WIGNER-DYSON UNIVERSALITY FOR SPECTRA OF REGULAR MATRIX 

PENCILS 

To study spectral properties of a regular matrix pencil H — AW we employ the Efetov super- 
symmetry approach p4|j35| |. A pedestrian introduction into the method can be found in pfl]. 

A convenient starting point is a representation of the spectral density p(A) of real eigenvalues 
A„ for the equivalent symmetric eigenvalue problem for the matrix H = W^^/^HW^^/^ in the 
following form: 

1 ^ 1 1 
p(\) ^ \^S(X~Xn)^T lim — - ImTr ^ (14) 

^^ ' N ^^ ^ ' ^c^+ottN (A±ie)l-H ^ ' 

which after a trivial manipulation can be rewritten as a derivative 

of a generation function Z±{J, A) defined as 

^±('^' ^) - det[(A±z.)W-H] ^^^) 

In a completely analogous way one expresses the two-point spectral correlation function ("cluster 
function", see Eq.(|3|)): 

1 

as 



J^2(Ai,A2) = (p(Ai)p(A2)) - (p(Ai))(p(A2)) = ;^Re/Cco„(Ai,A2) (17) 



1 



/Ccon(Ai,A2) = -rpr Imi lim 



a^ 



[(Z_(Ji,Ai)Z+(J2,A2))-(Z_(Ji,Ai))(Z+(J2,A2))] (18) 



N2 ,^+0.h_2~*0dJidJ2 

where the angular brackets stand for the ensemble averaging and we assume that Ai ^ A2. The 
expressions above are actually valid in the limit N —> 00, where one can show that 

32 



lim lim 



KZ±(Ji,Ai)Z±(J2,A2))-(Z±(Ji,Ai))(Z±(J2,A2))]=0 



e^+O ,/i.2^0 dJidJ2 

and we used this fact in Eg. (p7|) . 

To facilitate the ensemble averaging we represent the ratio of the two determinants in Eq.(|l9) 
as the Gaussian integral 



N 



N 



Z±(J,A) = (-1)^ / nM*.(±)]exp{±-^vl/t(±) W,,i\±ze + Jk)-H,, ^,{±)} (19) 



over 4-component supervectors ^i(±), 



,V.{±) 




,rf*,: 



drtdr* 



(20) 



with the components ri(+),ri(— ); i — l,2,...,iV being complex commuting variables and 
Xi{'^)iXi{^) forming the corresponding Grassmannian parts of the supervectors ^i(±). A 4 x 4 
diagonal supermatrix k = diag(0, 0, 1, 1) takes care of absence of the source J in the denomenator 
of the generating function Eqdlfl) . 

Since we are dealing finally with the averaged product of two such generating functions it is 
convenient to introduce the 8— component supervectors 



$, 




(21) 



and finally <& — ($1, ..., ^n) as well as the supermatrices A — diag(l, 1, 1, 1, —1, —1, —1, —1) and 
Jfc = diag(Jifc, J2fc)- 

It is also useful to remember that we expect to have non-trivial spectral correlations on the scale 
comparable with the mean spacing between neigbouring eigenvalues, i.e. when Ai — A2 oc 1/iV. 
Correspondingly, we introduce: Ai = A — s/2N , A2 = A + s/2N and consider s to be of the 
order of unity henceforth. All this allows us to write the product of two generating functions as 



^^■'(Jl,^2)=Z_(Ji,Ai)Z+(J2,A2) 



(22) 



d$exp<j — $'''A«)H$ + -$^ 



XA-{e-is/2N)l + Jk 



.W$ 



This expression can be easily averaged over the Gaussian distribution of H by the chain of 
identities: 



e 2 - 



&t / 



A®H*\ ^ ^-fj^StrA^ ^ 



dQexp 



-^StrQ2_i^$tAQ$ 



(23) 
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where a determines the variance of matrix elements of H, sec Eq.(O) and A = 
A^/^ (X^i *^j ® *^0 A^/^. The last relation trading the term in the exponent quartic with respect 
to $ for an auxilliary integration over the set of supermatrices Q is known as the Hubbard- 
Stratonovich transformation and plays a cornerstone role in the whole method. After substitut- 
ing Eq.(p3|) back into averaged Eq.(p2|) and changing the order of integrations one performs the 
(Gaussian) integral over $ explicitly. It turns out that in order to justify all these operations 
formally one has to restrict the supermatrices Q to a manifold paramatrized a,s Q — T~^PT, 
with P being a block-diagonal Hermitian supermatrix and T belonging to a certain graded coset 
space U0SP{2, 2/A)/UOSP{2/2) ® UOSP{2/2). A detailed discussion of this fact and an exphcit 
parametrisation of the matrices T could be found in [B4H3^] . 

The resulting expression for the averaged product of the generation functions turns out to be 
dependent only on eigenvalues Wi , i — I, ..., N oi the matrix W (this fact can be traced back to 
the "rotational invariance" of the Gaussian Orthogonal Ensemble (GOE) formed by H) and has 
the following form: 

dQexp --StrQ2_^Strln{(Ai-w.Q) + (ie + .s/2iV + Jfc)A} (24) 

where Vi = a/wi. 

Since we are interested in the limit (e, J, 1/A'^) — > 0, we can expand the logarithm in the exponent 
correspondingly : 



{{' 



Strln Al- 



+ {ie + s 



/2N + Jfe) a} 



Strln Al 



(Xi - v,Q\ + Str (Xl - v,Q\ a (ie + s/2N + Jk) - ^ Str f Al - v,Q\ 



-1 
v^Q ) A Jfc 



(25) 



so that differentiating the exponent over the sources Ji^2 yields the preexponcntial factors: 

N 



lim — — 

Jl,2^0 oJi 



exp{...} ex 2^ Str 



i=l 



"k^ 



and 



limji 2^0 



TV 



dJidJ: 



■exp{...} ex 2X^^=1 Str 



Al 



(Ai - w4) 



v^Q) '^^(ai 



1 



;a-i 



■ELStr 



Xl-v,Q] fcA±i 



^N 



e:Li str 



(ai - ihQ) k^ 



(26) 



(27) 



In the limit iV 3> 1 each sum over i is of the order of N, so that a contribution of second term in 
the expression above is by factor N larger than the contribution of the first one. In other words, 
we could restrict ourselves by terms linear in sources Ji,2 in the expansion Eq.(pq) above. We will 
make use of this fact later on when considering a more complicated situation in the next section. 
Taking all this into account, we arrive at the following integral representation: 
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1 9^ 

JC{Xi,X2) = -rp^ lim hm -— — — (Z_(Ji,Ai)Z+(J2,A2)) 



7V2 e^+0.h_2^QdJidJ2 

1 



4A^2 

AT 

Est 



lim e ^ / dQ exp 



AT 

X > Str 

1=1 



(Ai-«,Q) '^^^ 



NC{Q)~-{e-is/2N)^Str(\i-ViQ) A 

-i,A-i 



(28) 



Af 



EStr 



(ai - v.q) 



where 



£(Q) = ^ Str Q' + ^ E St^ In (ai - v,Q 



(29) 



In the same way one obtains the expression for the mean spectral density: 

N 



1 



p^^^-w^H^"^ h^QT.^'' 



-1 - A + 1 
Al— fiQ) k — - — 



exp 



^NCiQ)--eJ2Str[Xl-v,Qj A 

i 

(30) 



In the limit N —^ oo the integrals over Q in the expressions above are dominated by saddle 
points of the "action" C{Q) Eq.(E3) which satisfy the equation: 

Q = 4E--^^^- (31) 

relevant solutions of which belonging to the integration domains Q — T^^PT are parametrized as 
follows: 



,p ^tl + iqf-^kf 



where the real parameters t, q satisfy for q ^ Q the system of two equations: 



t = 



-E 



UijX - Vit) 

2 _i_ „2„2 



N ^ {\- v,tf + q 



and 1 — 



-E 



N ^^ {X- vay + q^v. 



(32) 



(33) 



Let us note, that the saddle-point solutions Eq. (|32|) form for q ^ a, continuos manifold 
parametrized by the supermatrices T. 

Using these expressions it is easy to invert the matrix I Al — ViQs.p) and check that 

N 



^Str 

i=l 



-1 - A± 1 
[XI- ViQs.p] k — - — 



(ai - ViQs.p] 



±2 



A - Vit 



*9l 



^StrfAl-WiQ^.p) A = iStr f-iATA ^ 



{X-v.,tf + q^vf ^{X-v.ty + qHf 
Viq 



Str 



f^^ATk 



1±A 



(A - Wit)2 + g2„2 

So the problem amounts to substituting these expressions into the integrand of Eqs.(28 30) and 
to performing the remaining integration over the coset space parametrized by the matrices T. The 
way of doing it described in much detail in |3J-|3^ and we refer the interested reader to those 
papers. 



(34) 
(35) 



12 



Here we mention only a few most important moments. First of all the integral in Eq.(BO|) is 
given by the so-called PSEW (Parisi-Sourlas-Efetov-Wegner) theorem due to a specific symmetry 
of the integrand ( in the limit e ^ the integrand contains only a part of the integration variables 
due to the projector (A + l)/2) The resulting mean eigenvalue density is merely given by: 



P(A) 
The system of equations Eqs(p3 






E- 



(36) 



(A - t;,t)2 + g2„2 

turns out to be in fact equivalent to a particular case of general 



results by Pastur and Girko |16|. 

Substituting the expression Eq. (|36|) for the mean spectral density to Eqs.(p4|) and the resulting 
formulae further to the integral Eq.(g8 29) for the averaged generating function one can express 
the "connected" part of the two-point correlation function /Ccori(Ai, A2) as: 



^con(Al,A2) — 



dn{T) Str 



T-'^ATk 



1 + A 



Str 



f-^kfk 



1-A 



exp 



iiTp{\) 



sStr 



(37) 

f-^kfk 



with dfi{T) being an appropriate invariant measure on the coset space. When deriving that ex- 
pression from Eqs.(28 29|j34| ) we noticed that the "disconnected" part of the correlation function 
is again given by PSEW theorem and cancels exactly the contribution from those few terms in the 
integrand which are not proportional to the mean density p(A). 

The expression Eq.(p7|) is our main result and is quite remarkable: its form exactly coincides 
with the corresponding expression for the two-point correlation function of random matrices from 
GOE, see e.g. [^. This fact was first demonstrated by Efetov who managed to perform a non- 
trivial integration over the coset-space explicitly ||34| and found that it reproduced the famous 



Dyson expression |3q] for the two-point function |28|: 

2 
3^2(Ai,A2) 



-t- 



sm s. 



dt 



sin (sgi) 



where Sg = T:p{X)s is a spectral distance s/N measured in units of the local mean spacing A = 
l/[NpiX)]. 

In the end of the present section we just mention that our method of deriving Eq.(p7[) can be 
easily generalized to matrix pencils whose random part H is an arbitrary matrix (real symmetric 



or Hermitian) with i.i.d. entries following the lines of papers [E2 26 1. In the same way one can 
show that all higher correlation functions also coincide with GOE expressions. 



III. SPECTRAL PROPERTIES OF A "FULL CONNECTIVITY" LC NETWORK 



Now we are well- prepared to calculate the mean density of resonances and the corresponding two- 
point spectral correlation function of the "full-connectivity" LC network as defined by eqs.(||,|lO|) 
Actually, it is convenient to rescale both matrices in Eq.|| as: H — > ^l/2 'ii, W 



1, w 
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This transformation obviously does not change generahzed eigenvalues of the pencil H — AW but 
facilitates bookkeeping of the leading terms in various expansions. 

In our consideration we follow the pattern of the previous section and introduce the generation 
function identical to Eg. (|2^) : 



i=l 



N 



I i<j 



where we used that 



for any supermatrix B as long as $]i?$j = $Ji?<i>i. We also envisaged that for the full-connectivity 
network the resonance density is nonvanishing as long as A ~ 1/N^''^ and introduced the scaled 
variable r — XN^''^ (see the Introduction). Then the typical spacing between neighbouring reso- 
nances should be of the order of A ^ N^'^/'^ . Since we can expect a nontrivial spectral correlation 
on the distances S* = A2 — Ai comparable with A we introduced the quantity e = e — is/{2N'^/'^) 
and consider s to be of the order of unity. All other notations coincide with those in Eq.([2^). 

Now we have to perform the averaging over h^u, (m? ^) = ^i -S, 1, ..., N , each taking values ±1 
with the equal probabilities. This is done in the limit iV — > 00 as follows: 






K'j 



n (1 - ^ [(*! - *]) ^ c^^ - '^^)] ' + ■■] - ^-pi-i^ E [(*! - 1']) ^ ^^^ - *.)] '} (40) 



1<'J •■ - I,] 

and in similar way we can average also over hA,^^ hB,i- Introducing now the notations: 
/C($a, $6) = ^Air ($t _ $t^ A ($, - $fc) + [($1 - $t^ A ($„ - $ 



/($) 



Ari/2 



$^ ( iei 



]^^ + '^0'^-i^^*'^*^' 



we can rewrite the averaged generating function as: 



i y i ij J 



(41) 
(42) 

(43) 



xexp^jT^TTlEf*!-* 



iel + Jfc 



($, - $,) 



!<J 
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For the case of Gaussian real symmetric matrices with independent entries H considered in the 
previous section, the key point allowing us to make progress was the exploitation of the Hubbard- 
Stratonovich decoupling; see Eq.(p3|). In a similar way, we proceed here with using a functional 
generalization of the Hubbard-Stratonovich identity suggested by us earlier p2| in the context of 
studies of the sparse random matrices: 



exp 



^ J2 IC{'^^,'^J) \ = J 2?(.9) exp I -^ I d$,d$,.g(a>,)C($„ <^,)g{<i>,) + '-f^ 9{'^^ 



(44) 
where the kernel C{^a, ^b) is , in a sense, the inverse of a (symmetric) kernel /C($a, ^b)- 

I d*/C($a, *)C('J', $6) = 5{^a,^b) (45) 

and (5($a,$b) plays the role of a 5— functional kernel in a space spanned by the functions (?(<&). 
Some hints towards understanding of the identities Eq. (H,E3) are given in the Appendix. 

With the help of these relations one easily brings the averaged generation function eq.(B3h to 
the form: 

{Z-^'{J^,.h))= jv{g)eM~NC{g) + 5Ci{g)} 

C{g) ^ ^ J d'!>ad^bg{'^a)C{<i>a,<i>b)9{<^b) -\n J d^ei^^^^+f^"^ 



SC,ig)= In 



iel + Jk 



(*.-*.)}" 



/ n,rf<i>. exp E.(t5(*.) + /(*.)) 



(46) 
(47) 

(48) 



As is clear from these expressions, the term Ci{g) is only a small correction to the main term 
in the exponent of the functional integral when N — > cx), J, e — > 0. Therefore, the functional 
integration over g can be performed by the saddle-point method (cf.Eqs.( p8| , pO| )). The saddle- 
point configuration gs{^) can be found by requiring the vanishing variation of the "action" C{g) 
and satisfies the following equation: 



ff(1>a) 



(49) 



When deriving Eq.(|49|) we have used Eq.(p5[). It is completely clear that in the limit N ^ oo one 
can safely disregard the term /($&), Eq.(B3) as being of the next order in 1/A^ in comparison with 
g{^b)- From now on wc just put f — everywhere. 

Given the form of the kernel Eq. (^l|) , we managed to guess the following explicit solution to the 
saddle-point equation: 



gs{<Pa) = 4(r - Gi)($t A$J + 4zG2($l$a) + *(*IA$,)2 



(50) 



Substituting such an Ansatz to Eq. (H9) one can verify by a direct calculation that it indeed satisfies 
the saddle-point equation provided the real coefficients Gi , G2 are solutions of the system of two 
conjugate equations: 
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G2 + iGi^ f duexp|^u(r-Gi+iG2)-y| (51) 

G2 - iGi = / duexp<^ --u(r- Gi - iG2) - — 

The following few identities might be helpful when performing a verification of such a fact. Suppose 
we have a 4— component supervector ^ like that defined in Eq.(Effl) and consider a function 
F($) = i^(^'l'v[/) vanishing on the bomidary of integration. Then it is easy to verify that: 

f d-9Fi-9)^F{0) , f d^bi^i^b)Fi^b)--0 (52) 

the first identity being just a particular case of the PSEW theorem mentioned in the previous 
section, see e.g. Pqj . 

For further analysis it is very important that a solution for the system Eq.(Bl|) exists for arbitrary 
— cx) < r < 00 such that G2(r) > 0. We will see later on that the mean density of resonances is 



merely given by p{r) — -G2(r). For \r\ ^ 1 one can easily infer from Eqs.(51) that G-zir) developes 
a gaussian tail: 

G2(r»l)«(7r/2)i/2exp{-^} 

It is necessary to mention that the equations Eq.(^l|) virtually coincide with those emerging in a 
study of density of eigenvalues of a transition matrix on a randomly diluted graph in a limit of large 
connectivity performed by Bray and Rodgers [p9l, see their Eqs.(20),(21) and fig.l in their paper. 
Indeed two problems have many common features, and the coincidence is hardly accidental. An 
indication of a relation between the problems comes from the structure of the matrix H, see eq.(^, 
which is actually very similar to the structure of the matrix considered by Bray and Rodgers. 
Moreover, the resonances we study are eigenvalues of the matrix H = W^^/'^HW^^/^, and it is 
easy to write down explicit expressions for the entries of H because of the simple structure of the 
matrix W, see eq.(g). One finds that H = -^H + 0(1/A^^/^) which provides a direct link to the 
paper |g9[ in the limit N ^ 00. 

The most important consequence of the existence of such a solution gs{^a), see Eq.(|50|), with 
G2 7^ is actually the simultaneous existence of a whole continous manifold of the saddle point 
solutions parametrized as: 

3T($a)-5s(T$a), with f^kf^k (53) 

Indeed, from the invariance property of the kernel Eq.(pd|): /C(<i>Q, <!>(,) = /C(r$a,T<I>f,) and from 
the (pseudo)unitarity of the matrices T: |Sdetr| = 1 it follows that gs{T^a) must be a solution 
to Eq.(|49|) together with any given solution gs{^a)- However, if it hadn't been for the condition 
G2 7^ all these solutions would trivially coincide: gri^) = g{^) for any T defined as above, i.e. 
the symmetry of the solution would coincide with the symmetry of the equation itself. 

In the actual case G2 7^ presence of the combination ^\^a which is not invariant with re- 
spect to a transformation <f>a -^ T^a makes the symmetry of the solution Eq.(BO) to be lower 
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than the symmetry of the equation Eq.(M9). This is an example of a very well-known effect of 
spontaneous breakdown of symmetry. The phenomenon makes the situation to be less trivial and 
generates the whole manifold of the saddle-point solutions. Different nontrivial solutions are actu- 
ally parametrized by the supermatrices T which are elements of the same pO| graded coset space 
UOSP(2,2/4)/UOSP{2/2) <g) UOSP{2/2) which aheady appeared in the previous section, see a 
discussion after Eq.(p3|). 

It is simple to satisfy oneself that C{gs) — 0, and the same obviosly holds for the whole manifold 
gx- Therefore, the only term that renders the expression for the generation function to be non- 
trivial is SCi{gT)- In the limit (e, J, 1/A^) ^ we can expand Eq.([48|) as: 



SC^igr) = In J 1 + -^ ^ / d^^d^, ($, - <i>,)^ (^^ii + j) ($, - $, 



)gi(9r(*.)+ST(*j)) I (54) 



2 
where we introduced a supcrmatrix: 



^7V3/2 j d$ $t (^ili + J^ $ eiST(*) ^ ^^3/2 g^^ 



W [ill + j] 



Waf3^ d<^^a<i>ie^OTi'S>) (55) 



and used that fact, that it is enough for our purposes to keep only terms linear with respect to the 
source matrix J (see the discussion after Eq.(p7|)). 

To determine the matrix elements of W it is convenient to use the equation: 

G2^i^a + ^Gl^iA<^a = J d<^b (^iA$b) ($Ia$„) e*^^^*') (56) 

which follows from the saddle-point equation Eq.(^) and Eq.(pO|). Now we make a transformation: 
^a.b -^ T^a.b and easily find that: 

G2$tf tf $^ + iGi$lA$a = J d$b ($lA$b) ($iA$„) ets-^*") (57) 

8 
= Yl *a,a*l^^A„aA;3/3i^aaW;3a 

Q,/3=l 

where Kaa = ~1 for a being a "fermionic" index and unity otherwise. Comparing both sides of 
that relation (i.e. the coefficients of the quadratic forms with respect to the components of the 
supervector ^a) yields a relation between W, the functions G12 and the matrices T: 

W ^G2f-^kfk + iGik (58) 

where we used f^f = kf-^kf. 

Using this fact we find that the averaged generating function is expressed as an integral over the 
saddle-point manifold parametrized by the matrices T: 

{Z''''{Ji,J2)) - /"d^(r)exp|i^G2(r)Str('f-iAfA) -f ^TV^/^ Str Jfc (Gaf-iAfA + iGiAU (59) 
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whith dfj.{T) being the invariant (Haar's) measure on the graded coset space. Here we assumed 

2 



an infinitesemal negative imaginary part to be included in s. Remembering that Jk = Jifc^i^ + 



J2fciy^ and making the correspondence: 

Giir) — > T7 / 7"; rrr^ — ^r^ and G^ir) — > tt / 7^; rr^T ^r^ (60) 

we see that the obtained expression hterally coincides H] with the averaged generation function 
derived in the previous section, that means with the generating function Eq. ( |24| ) after restricting the 
integration to the saddle-point manifold Eq.(^2|) and expanding with respect to source terms. Since 
that generating function underlied the expressions Eqs.(|37|) for the two-point spectral correlation 
function as well as the expression Eq.(|30|) for the mean spectral density, we immediately extract 
those quantities for our actual problem. As we already mentioned above, the mean resonance 
density is given by p(r) = —G2{r), and the functional form of the spectral correlation function is 
given by the same Wigner-Dyson expression as before, see Eq. (J3q) , provided one uses p{r) for the 
actual mean resonance density. 

Finally, it is necessary to mention that using the same method it is straightforward to demon- 
strate that all higher spectral correlation functions will be given by GOE expressions as well. 



IV. CONCLUSION 

In the present paper we introduced the full-connectivity model of a disordered reactance LC 
network and found that its spectral properties in the thermodynamic limit N —f oo could be 
efficiently investigated in the framework of a version of Efetov's supersymmetry method [p4|-|36| 
exploiting a generalized Hubbard-Stratonovich transformation introduced by us earlier [ p2| . We 
also studied spectral properties of regular pencils of random matrices with i.i.d. entries. In all cases 
we were able to derive the mean spectral density as well as to characterize fluctuation properties 
of spectra. The models studied turned out to be a faithful representatives of the Wigner-Dyson 
universality class. 

Actually, we hope that the present paper may provide a convenient background for the regular 
investigation of fluctuations in electric properties of disordered LC (more generally, RL — C) 
networks, such as the conductance Yab defined in Eq.(g), on-site potentials Vi, etc. Actually, it 
is possible to consider an LC— model of a "banded" type (i.e. of a large, but finite range, see 
| p3[ and [^ ) and derive |^^ the corresponding d-dimcnsional nonlinear cr— model which should 
adequately describe properties of the realistic LC— networks [g|,0| , including the effects of Anderson 
localization [Q These issues will be addressed in forthcoming publications |4^ . 
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APPENDIX A: A DIGRESSION ON FUNCTIONAL HUBBARD-STRATONOVICH 



TRANSFORM, EQ.(44) 



In this Appendix we present a heuristic demonstration of the validity of the relation Eq. (H3) for 
some kernels /C. We are not pretending that we are able to specify under what precise conditions 
our formal manipulations to be true: this interesting issue deserves to be studied seriously and 
goes much beyond our modest goals. Instead, our strategy will be a pragmatic one: to illuminate 
the origin of relations of that kind. 

Actually, the fact that $ is a supervector rather than a usual vector is of little importance for our 
consideration, so one may think about it as a usual vector (or even scalar). The way to incorporate 
anticommuting variables in the consideration sketched below is described in the Appendix A of 

ii- 

Let us deal then with a real symmetric integral kernel IC{^a, ^b) = !^{^b, ^a)- Then, generally 
we expect to exist a set of real eigenvalues Xi, and orthogonal and normalizable set of corresponding 
real- valued eigenfunctions e^{^): 

/d$b/C($„$b)e,($b) = A,e,($„) , /"d$e,($)e,,($)=5^, (Al) 

such that the kernel K, allows the following representation: 

/C($a, $6) = J2 ^-P-(*a)e.($b) (A2) 

where the summation goes over all v such that Ajy 7^ 0. Then we can write a formal chain of 
transformations: 

/=expJ-i Y. /c(i>.,<j>,)i =nc^p|-^E^-^''(*^)^-(*^)| =n°^p|-Y (Ee.($,)j 

Again, we do not specify conditions when it is allowed to interchange the summations, products 
and integrations but just assume simple-mindedly that the sequence of transformations presented 
above does not lead to divergent expressions. 



(A3) 
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A form of the expression Eq.(A4) suggests to consider it as a definition of a functional integral 



going over the space of functions: 

5(*)-E^i^'^-^-('^) ' ^(5)^n^ (A5) 

Finally, we introduce the kernel C($a, ^h) as: 

C(-^a, '^b) = Y. A;^e^(<i>a)e^(<i>fc) (A6) 

which obviously satisfies: 



where the ^g($a, *f b) plays the role of a (5— function in the space spanned by functions <?($), Eq.(A5): 



Moreover, using Eq.(A5) it is straightforward to verify that: 



which finally allows to rewrite Eqs.(AS,A4) as: 



e^P I - ^ E ^(*., *.) [ = / 'Dig) exp I -^ / d'i>6d$ag('i>a)C($a, '^b)g{'^b) -i^ai'^A 

(A8) 

which can be brought to the form Eq.(H4) by a trivial scaling transformation. 

Let us finally note, that our actual kernel Eq.(|4|) is of a separable nature and, moreover, has only 
few nonzero eigenvalues. Evidently, in that case the formal manipulations Eq. ( |A3| ) are expected 
to be most harmless. 
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